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Abstract. 

In general relativity, quasi-local energy-momentum expressions have been constructed 
from various formulae. However, Newtonian theory of gravity gives a well known and an 
unique quasi-local mass expression (surface integration). Since geometrical formulation of 
Newtonian gravity has been established in the covariant Newtonian space-time, it provides a 
covariant approximation from relativistic to Newtonian theories. By using this approximation, 
we calculate Komar integral, Brown- York quasi-local energy and Dougan-Mason quasi-local 
mass in the covariant Newtonian space-time. It turns out that Komar integral naturally gives 
the Newtonian quasi-local mass expression, however, further conditions (spherical symmetry) 
need to be made for Brown- York and Dougan-Mason expressions. 
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1. Introduction 

General Relativity (GR) is a theory for describing the dynamics of metric g and all other 
classical matter fields. Unlike energy and momentum densities of matter fields which are 
well-defined, gravitational energy cannot be localized because of the equivalence principle. 
Quasi-local idea is to define gravitational energy-momentum associated to a closed 2-surface 
(see the review article [8]). Relativistic expressions of quasi-local mass, energy, and angular 
momentum have been formulated from various approaches. However, Newtonian theory of 
gravity exists an unique and physically well-defined expression of quasi-local mass. Since 
GR will yield Newtonian theory of gravity in a certain limit, it is reasonable to believe that 
all of quasi-local expressions should return to the Newtonian quasi-local mass expression 
in the same limit. Newtonian and post-Newtonian approximations are asymptotic to general 
relativity by using the Newtonian sequence method [6], and hence Jeffryes [7] use it to verify 
the Newtonian limit of Penrose quasi-local mass. Unfortunately, it can return to the Newtonian 
mass and energy only by using transaction between 2-surface twistors and 3-surface twistors. 
Dougan-Mason mass [4] is another expression based on spinors and it has weak-field limit by 
transferring 2-surface integral into 3-volume integral. It is not clear to us that these quasi-local 
expressions with associated spin frames can directly return to the Newtonian quasi-local mass 
expression in the Newtonian space-time. 

The post-Newtonian theory provided a scheme to obtain an approximate solution 
of Einstein field equations and equations of motion based on two limits: slow motions 
and weakly gravitational field. Though Newtonian dynamics can be recovered in this 
scheme, it is not a covariant approach since coordinate condition (harmonic coordinates) and 
background Minkowski metric are required. However, Carter el al. [3] established covariant 
Newtonian multi-fluids and (magneto)-elastic solid models for neutron star crust by using 
geometric formulation of Newtonian gravity which has been largely developed by Kunzle [5]. 
Furthermore, if one try to probe the Neutron star, this requires a more compatible treatment 
between relativistic and Newtonian theories since the inner core may need to be described 
in terms of GR and it is sufficient to use the Newtonian theory for the outer crust. The 
covariant Newtonian space-time J\f is considered as a. flat Euclidean 3-space £ direct product 
a one-dimensional time manifold t and characterized by a degenerated (contravariant) flat- 
metric 7 and a one-form dt. It should be stressed that the metric-compatible and torsion-free 
connection V cannot be uniquely determined [5] by 7 and dt, so it is not necessary to yield a 
flat connection, i.e., Vx a ^x b = Vx 6 Vx a - One can introduce the Newton-Cartan connection 
D, instead of the flat connection, to obtain a geometrical formulation of Newtonian theory 
[3]. The detail construction will be reviewed in the following section. 

Our purpose in this paper is to examine Komar integral, Brown- York (BY) expression, 
and Dougan-Mason (DM) mass in the covariant Newtonian space-time. In [3], it yields a 
four-dimensional metric g = — (c 2 + 2$) dt ® dt + 77 in a pseudo-Riemannian spacetime M. 
and when speed of light c approaches to the limit c — > 00, one obtains {g/c 2 , g' 1 , V} — > 
{dt, 7, D}. Here, $ is considered as the Newtonian potential and rj denotes the (covariant) 
flat-metric on S. This provides a covariant transition between relativistic and covariant 



Quasi-local mass in the Newtonian space-time 



3 



Newtonian theories. We first use the metric g to calculate the three different quasi-local 
expressions and it turns out that Komar integral can exactly yield the Newtonian quasi-local 
expression, but BY and DM's energies expressions give zero. The result implies that the 
measure of these two quasi-local energy requires the information of non-flat metric on £ and 
leads us to ask a question: can a purely Newtonian type of star distort the geometry ofT,, i.e., 
generate pseudo-Riemannian curvature on it? It has been mentioned that g is sufficient to 
represent the relativistic corresponding metric of the covariant Newtonian theory. However, 
the post-Newtonian theory indicated that a Newtonian star should generate a spatial non-flat 
3-metric in the leading order 0(l/c 2 ). Thus, the Newtonian limit of GR will actually yield a 
non-flat 3-dimensional hypersurface. It leads us to generalize g by including a leading-order 
correction (a covariant 3-metric h) in i] which gives g = — (c 2 + 2$) dt <g> dt + (77 + h/c 2 ). 
It is easy to show that g can also return to the covariant Newtonian space-time {dt, 7, D} in 
the limit c — > 00. The three quasi-local expressions will be calculated by using g instead of g. 

The plan of this paper is as follows. In Section 2, we give a review of covariant 
Newtonian space-time and Newton-Cartan connection which is associated to the Milne group. 
In Section 3, we consider a specific form of non-flat 3-metric h which consists of two unknown 
functions in the 2+1 decomposition and are determined by Newtonian Poisson equations. 
The solution agrees with weak-field of Schwarzschild solution when one imposes spherical 
symmetry on it. In Section 4, we present a detail calculation of Komar, BY and DM's quasi- 
local energy by using g and also find the asymptotic Newtonian spin frame. It should be notice 
that the convention is chosen (H ) in Section 4.3 due to the spinor algebra usage. 

2. Milne Structure of Newtonian Space-time 

Since Einstein developed the special theory of relativity, Euclidean 3-space and time are 
merged into a four-dimensional manifold called Minkowski space-time which is characterized 
by a flat metric g M and Levi-Civita connection V. Later on, general relativity (a relativistic 
theory of gravity) was then generalized to four-dimensional curved spacetime. Though 
general relativity gives the Newtonian limit at a slow-motion approximation, it is a coordinate 
dependent approach, i.e., it requires to introduce the specific coordinates. To amend this 
unsatisfaction, the covariant description of Newtonian space-time provides a more consistent 
transition between relativistic and Newtonian theory of gravity. Furthermore, some problems 
of astrophysics require both relativistic and Newtonian treatments, e.g. inner core and outer 
curst of a neutron star, so to facilitate the transition can also improve solving these problems. 

Newtonian (Galilean) space-time is considered as a direct product of a flat Euclidean 
3-space £ and one dimensional Euclidean time line t, i.e., four dimensional fiber bundle. 
Each fiber S t contains a Euclidean (contravariant) flat metric 7 and one-form dt where 7 
is a degenerate metric and cannot upper or lower indices, i.e. no metric-dual definition of 
tensor fields. It has been found that conditions V7 = and Vdt = with vanishing torsion 
cannot give an unique connection V in the Newtonian space-time [5]. Nevertheless, V can 
be uniquely determined by introducing an ether frame e which is a unit timelike vector field 
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d£(e) = 1. By using 

Vat/j = 9 a ^-r^ = o, 

the torsion-free connection components can be obtained as 

r> = f^Vffi " ^(a^)7^ + ^^M^ - *(a V (2) 

where are the components of covariant 3-metric rj satisfying 

77(e, -) = 0, = 5/ - e%. (3) 

Greek indices (a, f3, ■ ■ ■ = 0, 1, 2, 3) denote the components of tensor fields with respect to 
coordinate basis {d^} and co-basis {Ax u } and round brackets indicate index symmetrization. 
Since 7 and dt are fiat metric, there should exist an Aristotelian coordinate system 
{t,X\X 2 ,X 3 } in which 7^ = 5^5^ + 62V + Ws 1 ', t„ = V> and also the 
corresponding ether frame = <5 M with rj^ = 5^5^ + 5^ 2 5 v 2 + 5 M 3 ^ 3 . Equation (2) 
then shows that r 7 Q/3 vanish in Aristotelian coordinates and thus V is aflat connection, i.e. 
V^V*, - Va,V a „ = 0. 

Besides coordinate transformations, there exist gauge transformations corresponding to 
the freedom of choosing e. It is easy to verify that V is invariant under Galilean gauge 
transformations 

e — > e = e + b (4) 

where 6 M are constants and satisfy dt(b) = 0. This corresponds to the invariance of 
Newtonian mechanics under Galilean transformations. To establish a covariant Newtonian 
theory of gravity in terms of geometry, Carter et al. [3] consider the Milne structure which is 
preserved by time-dependent Galilean gauge transformations (Milne gauge transformations), 
i.e. 6 M = in the Newtonian space-time. It follows that 7 /i "Va I ,b = 0. In the Milne 

structure of Newtonian space-time, V is not gauge invariant and has the following gauge 
transformations 

F M — F^ _ nV-f f (5\ 
vol 7 - 1 - ua - 1 - ua i* i>v v a- \~V 

Linear acceleration a is defined by a := Veb and a^(t) = dx M (a). It is worth to stress that 
the Milne gauge-dependent connection is still a flat connection. 

In the Newtonian theory of gravity, gravitational field g can be described in terms of the 
Newtonian potential 

g=-(7^V 9 »«9 M (6) 

and therefore d£(g) = 0. If one believes that linear acceleration a and g cannot 
be distinguished from any local Newtonian observer (a Newtonian type of equivalence 
principle), it then follows that Newtonian mechanics should be invariant under Milne gauge 
transformations. To obtain this invariant property, it is necessary to introduce a Milne gauge- 
invariant connection D (Newton-Cartan connection) 

^ 'vol ' ua fj^tuta- (7) 
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and g must have following gauge transformations 

g ^ g = g - a, $^$ = $ + « (8) 

to compensate transformations (5), where a M = 7 M7 Va 7 a. According to (1) and (7), one can 
verify that D is metric-compatible and torsion-free connection Ddt = D^y = 0. We then 
further to calculate the curvature tensor 

R(0„, d v )d a := D d D dv d a - D dv D d)i d a (9) 
which gives non-zero curvature 2-forms 

n = -^(V 9 „g)%"Adt (10) 
and Ricci 1 -forms 

p a = i d ^ a = -t a (y 9 ^y&t 1 (id 

where i dii denotes interior derivative. To complete the geometric formulation of Newtonian 
gravity, one introduces the field equation 

P a =47rGp%t a dt (12) 

which is equivalent to the Newtonian gravitational equation 

7^V 9m V^$ = AnGp (13) 

where G is the Newtonian gravitational constant. p M denotes mass flux 4-vector and mass 
density p = pH^. In terms of familiar vector calculus on S, (13) can be expressed as Poisson 
equation V • V$ = AnGp. 

It is convenient to define a 3-dimensional Hodge map # on E t with #1 = e 1 A e 2 A e 3 . 
{e 1 } are orthonormal co-frames associated with 3-metric r\ and Latin indices = 1,2,3. 

Without losing any generality, one can simply consider p = p(x l ) and then the covariant 
expression of Newtonian quasi-local mass yields 

M= I p a t n dV 



= [ P a t a 



, , d#d$ = -— <b #d$ 
4ttG ./ Et AttG J s=m 



1 



g*n„ y/\a\e ab dx a dx b = 7T77 <j> V$ ndS (14) 



8nG J s * " v 1 1 u " 4nG 
where x a (a, b = 2,3) are coordinates on the closed 2-surface S and n is a space-like 
(dt(n) = 0) unit outward normal to S. e a b denotes Levi-Civita e-symbol. It is worth to 
point out that (14) is valid for an arbitrary distribution of mass density p(x % ) on E 4 and Milne 
gauge invariant. 
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3. Relativistic Correspondence 

It has been approved that most of relativistic theories of gravity should regain the Newtonian 
theory of gravity in the limit c — > oo. In the post- Newtonian approximation, the leading 
order of Einstein field equations yields Newtonian Poisson equation in harmonic coordinates 
and it turns out that Newtonian potential $ appears in metric components g 00 = — c 2 — 2$ 
and gij = (1 — 2$/c 2 )<5jj. Though $ in g 00 is sufficient to represent the Newtonian limit, it 
indicates that a Newtonian type of star will actually generate a non-flat spatial metric. Since 
the covariant Newtonian theory has been established, it provides a covariant approach from a 
corresponding relativistic theory. Carter et al. [3] considered the following four dimensional 
non-degenerate metric 

g = n - c 2 dt <g> dt (15) 

with its inverse g' 1 = 7 — c~ 2 e®e in a pseudo-Riemannian manifold M, where c 2 = c 2 +2$. 
It can be shown that c — > 00 gives 

g' 1 -*!, c~ 2 g ^ -dt(g)dt, V -> D, (16) 

where Levi-Civita connection V yields the Newton-Cartan connection D. Though (15) is 
sufficient to represent the corresponding relativistic metric of covariant Newtonian theory, 
it seems inconsistent with the post-Newtonian approximation which gives a non-flat spatial 
metric. Furthermore, BY and DM's quasi-local energies give zero since 7 and 77 are flat 
3 -metric. It leads us to generalize g by including leading order correction on 77. This 
generalization will obviously preserve the covariant Newtonian theory in the limit c — > 00. 

Since the calculation of quasi-local quantities is associated with a closed 2-surface S, it 
is useful to apply 2+1 decompositions on 3-metric which gives 

rj = h ® n + a, 7 = n(g>n + 5" (17) 

with dt(n) = a(n, — ) = 0. h — g(n, —) denotes the metric dual of the space-like unit 
vector n and a corresponds to the metric on a family of closed 2-surface S n with their 
induced contravariant metric a = a a/3 d a dp. Thus, the generalization of g in terms of 
2+1 decompositions yields 

g — — c 2 dt®dt + n®h + a + (h®h + n®h) + a 
= - c 2 dt ® dt + f\ (18) 

where n and a are order of 0(l/c 2 ) and may be considered as small deformations of S n and 
its normal h. We further restrict deformations of S n as follows 

n = —n, a = — a (19) 

c c 

where / and q are functions of x^. The deformation can reduce to the weak- field limit of 
Schwarzschild metric in isotropic coordinates when f — q — GM/r. Furthermore, it can 
also yield post-Newtonian metric when / = q = — $. One may consider a more general 
deformation, however, (19) should be sufficient to represent a deformation generated by a 
Newtonian star. 
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It is convenient to introduce spherical coordinates {r, 9, 0} on E t and therefore g can be 
written as 

g = -d 2 dt <g> dt + (1 + ^)dr <g> dr + (1 + ^) r 2 (d# <g> d# + sin 2 d0 <g> d0) (20) 

cr cr 

where / = /(t, r, 6>, 0) and g = q(t, r, 9, 0). To solve Einstein equations 

^ ® e a = ^ (t - l -g T(X a , X a )^j , (21) 

we consider a Newtonian star which is characterized by the stress-energy tensor T = 
p c 2 dt ® dt and after some tedious but straightforward calculations, the leading-order Einstein 
equations yield Newtonian Poisson equation in {r, 9, 0} 

1 9 /2 S$, 1 d . . n d<$> , 1 <9 2 $ , ^ 

-7T r IT + ^-^7^ 8111 ^ + 2 • 2 Q = ^G P (22) 
H ar ar sin 6> a6> cw r 2 sin 9 ocp 2 

with / = q = — $, which agrees with post-Newtonian approximation. 

In the following section, we calculate Komar, BY, and DM's expressions in terms of (20) 

and then compare them to (14) in the limit c — > oo. For simplicity, we consider $, /, g are 

only functions of r, 0, in BY and DM's calculations. 

4. Quasi-local mass in the covariant Newtonian theory 

It is well known that (14) is valid for an arbitrary closed 2-surface, it may simplify our 
calculations by choosing a 2-sphere at r = constant on E t and then n = d r . Thus, (14) 
gives 

M = STg £ §7 ds = iTe i 57 r2s, "* dM * (23) 

From Section 4. 1 to Section 4.3, we present detail constructions and calculations of the Komar 
integral, BY and DM's quasi-local expressions. 

4.1. Komar integral 

The Komar integral is a covariant expression for quasi-local quantities based on Killing vector 
fields. In terms of a time-like Killing vector field £ with Hindoo <?(£,£) = —1, the Komar 
integral associated to quasi-local energy E K can be expressed as 

~4 



£ *=-87cr d « (24) 

where £ = g(£, — ) and* denotes 4-dimensional Hodge map with *1 = e°Ae 1 Ae 2 Ae 3 . Though 
no time-like Killing vector field exists in spacetime with metric g, it does asymptotically 
approach to a time-like Killing vector e/c in the limit c — > oo, i.e. 

1 / 2 9$ 1 \ 

lim - £ e £ = lim ^dt <g> dt + 0(— ) = (25) 

c^oo c c^oo \ COT & I 
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where £ e denotes the Lie derivative with respect to e. By using the asymptotical Killing 
vector field e/c, the Komar integral yields 

e c r 

M K := lim — = lim — — — (h *de 

c^oo C 2 c-^oo 8lT& Jg 

= -T^f/"x is <26) 

where e = g(e, — ). Equation (26) gives exactly the same covariant expression as (14) without 
choosing a specific closed 2-surface. 

4.2. Brown-York quasi-local energy 

The BY's quasi-local energy E BY obtained by using Hamilton-Jacobi analysis gives [2] 

E BY = ^—j>(k-k )dS (27) 
*J s 

where k := is the trace of extrinsic curvature k^ u of a closed 2-surface S and k 
denotes the reference values of quasi-local energy. A possible choice of k is by isometrically 
embedding S into a flat 3-dimensional slice of flat spacetime as a reference space. From (20), 
the unit normal h of a 2-sphere at t, r = constant is given by 

n = (1 + t)dr (28) 
<r 

and then the mean curvature k yields 

2 2 f 2 da 

k = -^Y4V^fi) a = — + -f- - (29) 
To obtain k , we start from a 3-dimensional reference space with a flat metric 

r] = dR ® dR + R 2 (d9 ®d9 + sin 2 9 d<p ® d0) (30) 
By performing the coordinate transformations r = R — rq/c 2 , one obtains 

fj = ( 1 + 2 ^ g M dr <g> dr + -^^(dr ® d9 + d9 ® dr) + ^-{dr <g> d0 + d0 <g> dr) 
\ c 2 or J c z 09 c z d(p 

+ ( y l + ^)r 2 ( y d9®d9 + sm 2 9dp(g)dp). (31) 
<r 

Therefore, calculating the mean curvature of S at r = constant in the reference space yields 

^ 2 ^ 2 <9(rg) 2 dq ^ 1 <9 2 g _^ cot 9 dq ^ 1 <9 2 g 
r c 2 r <9r c 2 9r c 2 r d9 2 c 2 r 89 c 2 r sin 2 9 d<p 2 
By substituting (29) and (32) into (27), we obtain 

c— >oo 

2f 2d(rq) 1 d dq 1 d 2 q\ AC 

r r ar rsm& <7tf w r sin a</>^ / 
Moreover, since Einstein field equation gives / = q = — $, M^y then becomes 

87rG J s \dr r sin # <9# (96 1 r sin 



lim 




c— >oo 


C 2 


1 




8ttG 


i( 
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By comparing (34) with (23), it turns out that only when $ = $(r), i.e., spherical 
symmetry, (34) can return to (23). Furthermore, weak- field limit of Schwartzschild solution 
$ = —GM/r yields correct result 

M B y = I \ dS = M. (35) 
4tt J s r 2 

4.3. Dougan-Mason quasi-local mass (H ) 

Dougan and Mason [4] used the holomorphic or anti-holomorphic conditions instead of 
the propagation equations of Ludvigsen-Vickers mass on a closed 2-surface, so that the 
Hamiltonian depends only on the boundary values of the spinors on the 2-surface. The energy- 
momentum can be defined as 

pAAL := J_j^ F ML dx * Adx d = ^£ F ML m « m P dS (36) 

by using Nester-Witten two-form F— = F^~dx a A dx^ and then its components are 

F^f := -i(zf V„z£ - Z§V a Z§). (37) 

The spin frame Zj = (Aa, p. a) which is normalized by e ab \aPb = 1 such that any solution 
of \ A can be written as A^ = Zj\a where \a is constant. The indices are lowered by the 
inverse e A b of e AB = \ A p B - p A \ B . 

Generally, e AB ZjZ§ may not be constant on the 2-surface S. Therefore, a natural 
definition of e— := e AB ZjZ§ demands the spinors satisfy the holomorphic (40) or anti- 
holomorphic (41) equations on S. So given a constant e—, we can define the mass by using 
the norm of the energy-momentum where 

! 

DM 

If we write N = e—, the mass can be calculated as (Bergqvist, [1]) 

m l M = j^(P^P^-P^P&). (39) 

Holomorphic condition is 

S\ A = 0, =>. 5A + <rAi = 0, 3Ai + p'A = 0. (40) 
Anti-holomorphic condition is 

S\ A = 0, =>■ S'Ai + a'X = 0, S'Ao + pA^O. (41) 
Then the quasi-local energy-momentum integral for a 2-sphere S can be written as 

I(r s ) = P—\aXa!_ = ^ + p\v) - A (9Ai' + p'\v)}dS. (42) 

We use the above equation and integrate by parts, then we obtain a very simple expression 

I(r s ) = ^ j{p\i\v + pA A '}d5 (43) 



m 2 DM := P^P^SabSa^ = P^Pax. (38) 
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where C = 1 for the holomorphic spinors and C = — 1 for the anti-holomorphic spinors. We 
will use the holomorphic condition in the later calculation. 

We first express the four dimensional non-degenerate metric (18) by using convention 
(H ). Then, in terms of null frame we have 

fiV = + ~~ ™ii m v ~ m fjnu, 9^ = + n^r — m^m" — m^m" (45) 
where the null tetrad can be written by ether frame and gauge invariant t M 



I, 


= o A o A > 






A A 1 
= 


— el _L r M 






= i A i A i 


= 2 (^ju ~~ r V)' 


n n 


= i A i A ' 


_ 2c 2 ' 


(46) 




= A l A i 


= V2^ X ^ ~~ ^m)' 


m il 


A A' 







The null frame satisfied the orthogonal relations Z M n M = l^m^m^ = —1. Also, the ether 
frame e M and can be expressed in terms of null frame where 

h = \t% + **), r, = ^-n„ (47) 

e" = c(y + n"), r " = ^-rA (48) 

and it satisfied the orthogonal relations e M t M = l,r'V /1 = —1 in the covariant Newtonian 
theory. Thus we chose e M = (1, 0, 0, 0) and = (1, 0, 0, 0) in (t, r, 0, </>) coordinate and use 
the metric 

ds 2 = ^2 _ (1 + / y dr 2 _ (r + 9 )2dd 2 _ (r + 5 )2 gin 2 g d(j) 2 (49) 

c c c 

with leading-order correction term /(r, 9, 0) and g(r, 0, 0). By choosing a corresponding null 
frame according to (46)-(48), we have 

I, = [5, -(1 + |), 0,0] P = [c- 1 ,(l + ^)~ 1 ,0,0], 



m 



m = |(1 + ^), 0,0] n" = -(!+/)-!, 0,0] 

* = [O,O,-^,^sin0] m" = [0,0,^ )7S ^] 



(50) 



where r' :— r + g/c 2 . The NP spin coefficients have the falloff 

/t = z/ = r = 7T = e = 7 = 0(— ), a — X — 0, a = j3 = p = fi = 0(1), 

<r 

with respect to c. 

The directional derivative along m can be expanded as 

2 

5 = m f V „ =6 +| + O(i) (51) 
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where 

o 1 i 2 n o 

S: = (d e -—dt), 5:= --5. (52) 

V2r sin0 r 

The key NP coefficients that we will use later are //, p, /3 and after some calculations yield 

2 2 2 

1 P 1 l it l o/5 i 

P = -r " ^2 + O(^), » = - Tr - £ + 0( ? ), /? =/? + O(-) 



r rcr c 



where 

i 2 i ^ n n f 9 % Vlcote 2 g%Jg 

^=2' P-=9r9-f--, P--=^— > P := ~r P + &- ( } 

The asymptotic constant spinor A^ = Aio^ — A <vi can be expressed as 

2 2 

Ao=Ao+^ + 0(i), Al 4i+^ + 0(i). (54) 
By substituting the relevant asymptotic expansion terms into the surface integration, we 

have 

pAiAi' — /iA A ' 
ooo ooo 

= PAiAi' - ^AoAo' (55) 

1 20 20 02 02 00 2 00 2 1 

+ -^[PAiAi' - ^AoAo' + PXiXv - PXoXo' + PAiAr - /^AoAc] + 0(— ). (56) 

We will show that only (56) will survive. 

The first order 0(1) holomorphic conditions are 

00 00 00 00 00 

5Xo - PXo= 0, SXi + (3Xi - VXo= 0. (57) 



Therefore, the solutions are 



o o (A o ^ o o (hod) 

Ao= C_i A-i exp(-i£) + d Ai exp(i J), Ai=B-i exp(-i£)+ Pi exp(i£) (58) 



where 



A-i = sin-, Ai=cos-, (59) 

fl • 

o cos £ 4i o sin £ A i 

B - =c ^r = Cl vr Bl= - c -v# = - c -w (60) 

We substitute the first order spinors and the related NP coefficients into (55) and integrate 



it. We then find P 00 vanish. Furthermore, we can find P 00 =P n from symmetry of the 

oo 

solution and P 10 =P 01 = from f Q * exp(i0)d0 = 0. The second order 0(l/c 2 ) holomorphic 
differential equations are: 

o 02 n o 

(5 - (3) Ao -f Ao= 0, (61) 
2r 

o 2 02 n o Q q - f o 

(5 + P)\i- PXo +7T Ai + Ao= 0, (62) 



Quasi-local mass in the Newtonian space-time 12 
where 

2 2 02 22 02 

Ao=A-i exp(-i-)+ Ai exp(i-), Ai=S-i exp(-i-)+ Si exp(i-), (63) 
and 



2 r 1 " 2 f 6 

A-i =C-! 5gd9'A-i, Ai= Ci \ 5 gd9 A 

J 6 JO 



(64) 



2 1 f n H 5 qd9' A-i oo oo 
B-i = - -—^ / [C-i Jo y -+ 5 g B-i +C^(d r g - f) A-!] Ai M, (65) 



V2 Ai 



2 l [\ f n e 5 gd9' Ai o o , oo 

V2A-i Jo V2 

2 2 2 2 

where A-i, Ai, B-i, Bi at 0, ir vanish. 

02 02 0202 

We can further find § A?A^ dS = § \\\% dS = and X\X^,—X^X^,. Therefore, 
only (00') component will contribute to the energy-momentum surface integration. Thus we 
substitute the second order spinors and NP coefficients into the (00') component of (56), we 
then find 

P°°' = ^ j\d r g - f) sm(29) j\d r g ~ f) sm9 sm9'd9d9' 

and 

|N| = |AgAi-A?A5| = -^ + 0(^). (67) 
The mass expression (38) can be expanded as 

km m DM = „ (68) 

| N | 

= ic So ^ 9 ~ f) Sin(2 ^ 69 ~ 2^ So Se> ^ 9 ~ f) 6 ^' ^ & 
For Einstein field equation satisfied, we have / — g/r — — $(r, 9, (ft). If we consider spherical 
symmetry of g, f and / = GM/r = — $(r) and g = — r$(r), then the DM mass becomes the 
quasi-local Newtonian mass 

lim m DM = <f d r $(r) dS. (69) 

c^oo 47TG J 5 



5. Conclusion 



In this work, we verify that Komar integral can yield the Newtonian quasi-local mass 
expression without choosing a specific 2-sphere or referring to spherical symmetry of the 
Newtonian potential $, however, the Brown- York expression and Dougan-Mason mass can 
give Newtonian expression only in the spherical symmetry of $. Since general relativity has a 
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well known Newtonian limit, it is reasonable for us to argue that it should yield the Newtonian 
expression for all of quasi-local expressions in the covariant Newtonian space-time. We 
discover that the differences between Newtonian expression and Mby are contributed from 
reference of energy surface density k which is obtained by isometrically embedding a 2- 
sphere in a flat reference 3-space. Moreover, we find that the Newtonian source will distort 
not only NP spin coefficients but also the spin frame for the leading-order term of the 
Dougan-Mason mass expression. 

The Newtonian theory is sufficient enough to describe the most of the astrophysical 
problems. It would be an interesting question that can we measure the correct amount of mass 
for an arbitrary Newtonian source by using these quasi-local expressions? Unfortunately, 
we find that only for the spherical symmetry source Brown- York and Dougan-Mason mass 
can get the surface integration of the Newtonian mass. Whether this problem is due to the 
expressions themselves or some other technical problems requires a further investigation. 
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